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1. INTRODUCTION. 


THE classical method of treating the thermodynamic properties of 
gas mixtures is to combine from the beginning of the discussion the 
methods of thermodynamics and statistical mechanics. In _ this 
manner van der Waals,' and van Laar and Lorenz? gave relations 
for the properties of gas mixtures based on the ideas from which 
van der Waals’ equation of state was derived; and Keyes*® obtained 
relations by use of the Keyes‘ equation of state. 
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Recently there has been proposed another method of approach, 
which involves the separation of the thermodynamic and statistical 
treatments until the end of the discussion when the general equations 
derived by thermodynamic methods are integrated by use of an equa- 
tion of state for gas mixtures. Some assumption or assumptions are 
necessary as a starting point for the thermodynamic discussion and 
Gillespie® has shown that general relations for the isothermal variation 
of certain properties of gas mixtures can be derived from the assump- 
tion of Dalton’s law of partial pressures at low pressures. De Donder,® 
Van Lerberghe’ and one of the authors® have also treated this problem. 

It has been found® that two assumptions are required in order to 
derive expressions for the variation of the thermodynamic properties 
of gas mixtures with both temperature and pressure. The equations 
so obtained may be called the general thermodynamic relations since 
they depend only on the two assumptions and the general laws of 
thermodynamics. In order to evaluate the partial derivatives which 
each of these equations contains, it is necessary to have available 
an equation of state for gas mixtures and such a relation must neces- 
sarily depend to a greater or less extent on the methods of statistical 
mechanics. 

In the present paper an equation of state for gas mixtures based 
on the Beattie-Bridgemen equation of state is considered from the 
standpoint of the recent experimental data. 


2. Previous INVESTIGATIONS. 


We are concerned here with the methods that have been proposed 
for calculating, from data on the pure component gases, the pressure 
(or volume) of a mixture of gases when the volume (or pressure), 
temperature and composition are given. 


(a) Simple Additive Relationships. 


There are four additive rules for the treatment of gas mixtures 
which should be mentioned: 

1. Dalton’s law.'° This is usually stated: The pressure of a mixture 
of gases is equal to the sum of the pressures of the component gases 
existing each by itself at the temperature and in the total volume of 
the mixture. This may be called the law of the additivity of pressures 
at constant temperature and volume. It is well known that the law 
does not hold for the successive additions of portions of a pure gas 
to a given volume unless the gas follows Boyles’ law. 
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On the basis of Dalton’s law, the equation of state of a gas mixture 
would be obtained by adding the equations of state for the component 
gases in the following manner: 


= filV, T, ny) + f(V, T, Ne) + warn (1) 
In this relation p is the total pressure of the gas mixture; V, 7’, ni, 
no, -*: are the total volume and the temperature of the mixture and 


the number of moles of the Gases 1, 2, --- composing the mixture; 
and 
pi =Sfl(V, T, nj) = 1,2, - 


are the equations of state of the pure gases. 
When applied to ideal gases, Dalton’s law gives the result that: 


Pi = Pri 


where p is the total pressure of the mixture and 2; the mole fraction 
of Gasiin the mixture. The product pz; is usually called the “‘ partial 
pressure.’ | 

2. Gibbs-Dalton law. From a consideration of vapor pressures, 
Gibbs" proposed a modification of Dalton’s law: The pressure of a 
mixture of gases is equal to the sum of the pressures of the component 
gases existing each by itself at the temperature and with the same 
value of its chemical potential as in the mixture. Lurie and Gillespie” 
were the first to make a study of the Gibbs-Dalton law. They pointed 
out that this law leads to the same expression as Dalton’s law, equa- 
tion (1), for the equation of state of a gas mixture. 

3. Amagat’s law. From a study of his measurements of the 
compressibilities of oxygen, nitrogen and air, Amagat® concluded 
that: The volume of a mixture of gases is equal to the sum of the 
volumes of the component gases existing each by itself at the tempera- 
ture and under the total pressure of the mixture. This may be 
called the law of the additivity of volumes at constant temperature 
and pressure. Since volume is an extensive property, Amagat’s law 
evidently holds for the “mixing” of various portions of a real gas 
each being at the same temperature and pressure. 

Amagat’s law leads to an equation of state for a gas mixture obtained 
by adding the equations of state of the component gases in the fol- 
lowing manner: 


V = F,(p, T, m1) + Fo(p, T, ne) + --- (2) 
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In this relation V is the total volume of the gas mixture; p, 7’, m1, 




























na, --* are the total pressure and the temperature of the mixture 
and the number of moles of the Gases 1, 2, --- composing the mix- 
ture; and 


v= F (p, T, n;) = l, 2,° > % 


are the equations of state of the pure gases. 

4. Lewis and Randall rule. On the assumption that gas mixtures 
could be considered as perfect solutions, Lewis and Randall" proposed: 
The fugacity of a gas in a mixture is equal to the product of the 
fugacity of the pure gas existing at the temperature and under the 
total pressure of the mixture multiplied by the mole fraction of that | 
gas in the mixture. Gillespie’ has shown that this rule applies only 
when Amagat’s law holds for the gas mixture. Hence the Lewis and 
Randall rule leads to the equation of state (2) for gas mixtures. 

Thus the first two laws result in an equation of state of the type 
(1) for a mixture of gases, while the latter two give an equation of 
the type (2). Either Dalton’s or Amagat’s law can be used wtihout 
the aid of an equation of state provided that compressibility data 
on the pure gases are available. ' 

It has been generally accepted’® that Amagat’s law gives more . 
satisfactory results than Dalton’s law for the permanent gases. Lurie 
and Gillespie’? have shown that the Lewis and Randall rule gives 
more satisfactory results than the Gibbs-Dalton law when applied 
to the calculation of the “equilibrium pressure” of ammonia in 
mixtures of nitrogen and ammonia. Keves and Burks" found that 
Amagat’s law is a better approximation than Dalton’s law of “ partial 
pressures” in terms of weight fraction (rather than mole fractions) for 
mixtures of nitrogen and methane. Masson and Dolley” found ‘ 
that Amagat’s law gave a better representation of their measurements 
on argon-oxygen mixtures than Dalton’s law although the reverse 
was the case for argon-ethylene and oxygen-ethylene mixtures. All 
of these comparisons show that Dalton’s and Amagat’s laws each 
give an increasingly good representation of the experimental data as 
the pressure is diminished, but that they may lead to large errors at 
moderately high pressures especially when one of the gases is much 
more compressible than the other. 


(b) The van der Waals Equation of State. 


There have been many investigations concerned with the repre- 
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sentation of the pressures of gas mixtures by use of van der Waals’ 
equation of state:'® 
 —_— 


v—b Yr 


P= (3) 
in which R is the ideal gas constant, a and b constants for each gas, 
and v the volume of a unit mass, say one mole, of the gas. 

General considerations indicate that the equation of state of a 
mixture of gases which do not react chemically should be of the 
same form as that for the pure component gases and the same con- 
clusion is reached by use of statistical mechanics.'? Thus from the 
standpoint of empirically representing the compressibility data, it 
should be possible to select values of a and 6 which give as satis- 
factory a result for a mixture of gases as for the pure constituent gases. 
Since the time of van der Waals it has further been the expectation 
that the values of the parameters for any given gas mixture could be 
calculated from the constants of the pure component gases and the 
composition of the mixture. General considerations and statistical 
mechanics again concur in indicating that each of the constants a 
and 6b should be combined with the composition in the most general 
quadratic manner. Thus for the b constant we may write: 


b, = bay? + Qhiox 102 + bow 2" (4) 


In this relation b; and /». are the values of the 6 constant for the pure 
Gases 1 and 2, b, the value of this constant for a binary mixture 
of these gases the composition of which is expressed in terms of the 
mole fractions 2; and 22, and by» a so-called “interaction constant” 
which arises from the interaction of unlike molecules. Thus when 
applying equation (3) to the gas mixture we can replace b by 6, and 
a by a,, using for v the volume of an average mole of mixture. 
Several methods for evaluating the interaction constants have been 
suggested, the two simplest of which are that they may be taken 
as the arithmetic, or as the geometric mean of the corresponding 
constants for the pure gases. Thus for the b constants we may write: 


bye = Ya(hi + re) (5) 
or: 
big = (bybe)'” (6) 


When these are substituted into (4) there results for 6,: 


h. = bya, + bere (7) 
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61/2 —_— 5/22, a bo! /*29 (8) 


The first of these methods for computing b,, equation (7), will be 
called linear combination of constants and was used to some extent 
by van der Waals and by Berthelot”® for the 5 constant; the latter 
method, equation (8), will be called linear combination of b' or linear 
square root combination and was used by Galitzine”! and by Berthe- 
lot®° for the cohesive constant a. 

On the assumption that the molecules are elastic spheres and that 
the distance of closest approach of the centers of any two molecules 
is the sum of the radii of the two molecules in question, Lorentz” 
derived for the bh, interaction constant the value: 


byt! = 16 (b¥3 + dQ") (9) 
Thus: 
b, = byxy? + 14 (O13 + be!?)3 x 22 + bore? (10) 


We shall call this the Lorentz method of combination. It has been 
quite generally employed* for b,. It is questionable however if the 
b constant should be linked too intimately with the distance of 
closest approach of the molecules. On the basis of Phillips’™ view, 
any kind of interaction between the molecules of a gas would give 
rise to a b term; that is, the pressure of the gas would be increased 
by any type of interaction, regardless of whether the interaction is 
due to attractive or repulsive forces, regardless of the law of the forces 
and further, independent of the size of the molecules, which may be 
infinitesimal. 

Because of the failure of van der Waals’ equation to give even an 
approximate representation of the compressibility of pure gases over 
a wide range of temperature and density, it has never been applied 
to compressibility data on gas mixtures except” at relatively low 
pressures. 

(c) The Keyes Equation of State. 


The first quantitative treatment of gas mixtures at relatively high 
pressures was by means of the Keyes equation of state:* 
RT A 
v— Be~*"—— (n+ 1)? 


in which R# is the gas constant and a, 8, A and / constants for each 
gas. Lurie and Gillespie’® applied this equation to the calculation 


(11) 
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of the equilibrium pressures of ammonia in a mixture of ammonia 
and nitrogen to 60 atmospheres at 45° C. They used linear combina- 
tion for each of the constants R, a, 8, A‘? and 1 as suggested by 
Keyes and stated that the computed values agreed with the experi- 
mental data possibly to within the experimental data. Keyes and 
Burks’ measured the compressibility of mixtures of nitrogen and 
methane to about 300 atmospheres and from 0 to 200° C. They 
determined the values of the four adjustable constants of (11) for 
each of the mixtures investigated, thus definitely proving that the 
compressibility of the mixtures and of the pure gases can be repre- 
sented by the same type of equation of state. They further found 
that each of the constants as determined from the data on the mixtures 
could be satisfactorily calculated from the values of the constants 
for the pure gases and the composition by the method mentioned 
above. 
(d) The Beattie-Bridgeman Equation of State. 


One of the authors® has applied to gas mixtures the Beattie- 
Bridgeman” equation of state: 





p= ait= e+ B)— = (12) 
A = Ao(1 — a/v) 
B = B,(1 — b/v) 
e = c/vT? 


in which R is the ideal gas constant and Ag, a, Bo, b and c are constants 
for each gas. The pressures of the nitrogen-methane mixtures as 
measured by Keyes and Burks are well represented when the values 
of the parameters for the mixtures are computed by linear combination 
of each of the constants Rh, Ao,'/? a, Bo, b and c. The values of the 
constants of the pure gases used were those determined from the 
data of Smith and Taylor” on nitrogen and of Keyes and Burks”® 
on methane; the former having used the same apparatus for the 
compressibility measurements on nitrogen as the latter used for 
methane and the mixtures. 


(e) Other Methods. 


Several other methods of attack on the problem of the equation of 
state of gas mixtures have been made. Thus Lennard-Jones and 
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Cook** have shown that the second virial coefficient should be quad- 
ratically combined in the manner of equation (4) but that the inter- 
action constant is not related in any simple manner to the coefficients 
for the pure gases. A relation for the interaction constant is given 
in terms of the force constants appearing in the equation representing 
the attractive and repulsive forces between molecules of unlike type. 
The relations obtained are too involved for general use. 

Deming and Shupe,*® extending a method suggested by Kleeman," 
have calculated the pressures of the binary mixtures of nitrogen and 
methane investigated by Keyes and Burks by combining each of 
the coefficients in the general virial equation: 


Ay Ao Az Ay 
PS Ss a. (13) 
v vy v v 
according to the method: 
Au: = (Dasemi(EmsAmi)!/™]” (14) 
where 4A,,; is the mth virial coefficient of the mixture; Am;, Ams, -*: 
the mth coefficients of the pure Gases 1, 2, ---; 21, 22, --- the 


mole fractions of the gases in the mixture; and ¢«,; = +1 or 
— 1, whichever is required to make the product ¢,;4m; positive. 
The summation extends over all the gases in the mixture. The 
method is applicable when the values of the virial coefficients for 
the pure components have been determined at a single tempera- 
ture only. Deming and Shupe found that this method does not give 
quite so good an overall representation of the nitrogen-methane data 
as the Beattie-Bridgeman equation of state for mixtures although in 
the low temperature high pressure corner of the table the representa- 
tion is somewhat better. 

Van Lerberghe*® has noticed that the generalized equation of state 
for a gas mixture: 


p=f(T,V, mi, ne, +++) (15) 


in which V is the total volume and m, me, --- the mole numbers of the 
Gases 1, 2, --- in the mixture, can be developed into a Maclaurin 
series. This gives for a binary mixture the result: 


iy 
p= pt pt ron (16) 


in which p; and pe are the pressures exerted by the Gases | and 2 
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respectively at the temperature 7’ of the mixture and the molal 
concentrations ;/V and n/V. The quantity wy is a function of 
the temperature and must be determined for each pair of substances. 
This relation represents the methane-nitrogen data well but it should 
be noted that in addition to data on the pure gases it is necessary to 
have available information concerning w 2: at each temperature, or 
the empirical temperature function for wy for each pair of gases. 
Thus new constants which cannot be calculated from those of the 
pure gases are introduced. This limits the usefulness of the method. 


3. Previous WorK ON THE METHODS OF COMBINATION OF THE 
CONSTANTS. 


Three methods of combining constants have been mentioned: 
linear, equation (7); linear square root, equation (8); and Lorentz, 
equation (10). For any given values of the constants for the pure 
gases, linear combination leads to the greatest value of the corres- 
ponding parameter for the mixture, linear square root to the smallest, 
and Lorentz combination to an intermediate value. 

When an equation of state constant has the dimensions of volume, 
linear combination for the corresponding p#iameter for the mixture 
gives the simplest resulting equation; when the constant has the 
dimensions of volume squared, linear square root combination gives 
the simplest result; while the Lorentz method does not lead to simple 
equations in any case. In the first two cases mentioned the constants 
can be combined with the mole-numbers of the component gases 
instead of the mole fractions. Many of the uses to which an equation 
of state for gas mixtures Is put require differentiation with respect 
to mole-numbers, and the resulting equations are simpler when 
mole-numbers can be used in the combination equations (7) or (8). 

There have been three methods proposed for determining the 
interaction constant b,. and two for aj, the Lorentz derivation not 
applying to this latter case. It is for experiment to determine: 

(1) Which of these methods of combining constants gives the 
most satisfactory representation of the experimental compressibility 
data on gas mixtures. 

(2) Which of these methods, if any, leads to an equation of state 
for gas mixtures that, when used to integrate the general equations 
for the thermodynamic properties of gas mixtures, gives relations 
which represent these properties satisfactorily. 

(3) Which of the three methods: combination of the constants in 
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an equation of state, Dalton’s law, Amagat’s law gives the best results; 
and how much better is the preferred method. 

The present paper deals with the first and last of these statements. 
There have been several investigations that have treated the question 
of what is the best method of combination for a given constant. 
We are concerned mainly with the a and 6 constants of van der 
Waals’ equations or the corresponding constants A and £8 of the 
Keyes or Ag and By of the Beattie-Bridgeman equation. The methods 
of combination used for the higher order constants in the latter two 
equations of state are not so important since in general these constants 
do not have so great an effect on the calculation of the pressure. 

In the evaluation of some of the thermodynamic properties of gas 
mixtures, the cohesive constant a has a relatively greater effect than 
any of the others. Eucken and Bresler®* from their own measure- 
ments and those of Pollitzer and Strebel* on the effect of an indifferent 
gas on the vapor pressure of liquids calculated values of the inter- 
action constant aj. by the use of certain expressions for the equilibrium 
pressures and the van der Waals equation of state. The /. inter- 
action constant was calculated by the Lorentz method, but it was 
remarked that the method of combination used for 6 did not greatly 
influence the result since the 6 constant does not enter markedly 
into the calculation. They found that the best representation of the 
interaction constant was given by linear square root combination 
for a. 

Gillespie* investigated the limiting value of the expansion of a gas 
in a gas mixture as the pressure approaches zero at constant tempera- 
ture and composition. The general relation for the intercept is: 


Ao: — 2Aor2 + Aw 
RT 


where ?, is the partial molal volume of the Gas | in the mixture and 
v, the molal volume of Gas | in the pure state at the temperature and 
under the pressure of the mixture. The same expression results 
whether the van der Waals, Keyes, or Beattie-Bridgeman equation 
of state is used provided that in the latter two equations all of the 
higher order constants are combined according to the linear method. 
Since Gillespie actually used the values of the constants of the latter 
equation of state in his computations, it is written here in terms of 
these constants. 

Twelve binary mixtures for which compressibility data existed in 





Limit (7; — 7) = 2¢? 
p—od0 


— (Bo. — 2Bor+ Ba) | (17) 
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the literature were investigated by Gillespie; five for argon-ethylene 

mixtures, three for oxygen-ethvlene, one for hydrogen-nitrogen, and 

three for nitrogen-methane. He found that the value of the intercept 

obtained by extrapolation of the values of v;— v; determined from the 

experimental data to zero pressures is not zero as had been assumed 

by some writers but has in all of the cases studied a finite positive 
value. 

When linear combination, equation (7), is used for Bo, the quantity 
(Bo: — 2 Boe + Bee) is zero [by equation (5)] and when either linear 
square root or Lorentz combination, equations (8) or (10) respectively, 
are used this quantity is positive since twice the interaction 
constant Boy determined by (6) or (9) is less than the sum of the 
constants for the pure gases. Thus linear combination for Ao is 
excluded since in this case the contribution from the A» constant to 
the value of the intercept would be zero and the calculated intercept 
would be zero or negative. Moreover Gillespie showed that when 
linear square root combination was used for Ao, and either linear or 
Lorentz combination for Boz, a good representation of the experi- 
mental intercepts was obtained. The argon-ethylene data pointed 
to the Lorentz method for Bo,, the oxygen-ethylene to the linear 
method, while the other two binary combinations did not strongly 
support either method but perhaps favored the linear method slightly. 
In no case was it impossible to exclude either method since the effect 
of Bo on the calculation was smaller than that of the Apo constant. 
On the other hand the linear square root combination for Bo, seems 
to be excluded since in all cases in which any doubt existed the values 
of intercepts obtained by extrapolation of the plots lay rather between 
those calculated by linear and Lorentz combination for Bo, whereas 
linear square root combination for Bo, would cause the calculated 
intercept to fall below the Lorentz value. 

This investigation definitely indicates the linear square root combi- 
nation for Ap and either linear or Lorentz combination for Bp with 
the linear method slightly favored. 


4. THe EQUATION OF STATE FOR MIXTURES USED IN THE PRESENT 
COMPARISONS. 


The virial form of the equation of state (12) when written for a 
gas mixture is: 


af (18) 








-™ 
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Bn = R,, TBon — Aon os R,c,/ T2 
Yn ot R,, T Bony, + AonGn sm R,Bonn/ T? 
bn _= R,Bonbran 7° 

















where R,, Aon, Gn, Bon, b, and c, are the values of the parameters for 
a mixture consisting of m1, m2, . . . moles of the Gases 1, 2, . . . respec- 
tively and V is the volume of =n; moles of the mixture. Any equation 
may be used for computing the values of the parameters from the 
constants of the pure gases and the composition; but a method that 
gave good results in the calcilation of the pressures of nitrogen- 
methane mixtures* and of the expansion*®® of gases on mixing at low 
pressures is linear combination for each of the constants R, Ao'”, 
a. Bo, b and c, namely: 






























R,=nRkR+nR+--- = Di(n;)R (19) 
Aon = (my Agi? + rede’? + +--+)? = (ZniAoi‘*)? 
An = NA, + Note + --- = >(n,a;) 
Bon = m1 Bo, + neBr + --- = =(n:Bo:) 
b, = nyby + nebo t+ --- = D(n;b;) 
Cy = NC, + Nolte + °°: = L(n;) | 
In these relations 1), #2, -:- may represent: (1) the number of moles of 


the component gases, in which case the volume V used in (18) is the 
total volume of Yn; moles; (2) the mole fraction of the component 
gases, the volume v being the volume of an average mole of mixture, 
that is V/n;; (3) the weight fractions of the component gases, the 
volume v referring to one gram of mixture. In the latter case R has 
a different value for each gas and R, = =(n;h,j). 

The method of combination (19) not only gives a satisfactory 
representation of the experimental data in the several cases in which 
it has been used, but leads to relatively simple equations for the 
thermodynamic properties of gas mixtures. Provisionally, we may 
consider that the method of combination is settled for all of the 
constants except for By which, in addition to the method shown in 
(19), may also be combined either by the Lorentz or the linear square 
root method. In the former case all of the equations of (19) remain 
the same except that for Bo, which becomes: 


Bor = rP Bou + 4rire(Ba'? + Bo)? + 22°Be (20) 
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and in the latter case it is necessary to replace the relation for Bo, 
of (19) by: 
Boz = (22;Bo;''")? (21) 


In (20) and (21) only the mole fraction may be used since the quantity 
By has the dimensions of volume. For the calculation of the pressures 
of a gaseous mixture, the use of either of these methods of combination 
does not introduce serious difficulties; but either expression leads to 
rather involved relations for some of the thermodynamic properties 
of gas mixtures, especially those requiring differentiation with respect 
to the composition. 


5. THe EXPERIMENTAL DatTa. 


The earlier work** on the compressibility of gaseous mixtures was 
rather fragmentary. Recently, however, there have been several 
investigations that extend to moderately high pressures: 

1. Kuenen, Verschoyle and van Urk*’ measured the compressibility 
of two mixtures of oxygen and nitrogen. ‘These data are not treated 
in the present paper since they are for the most part in the critical 
region and also because the values of: the equation of state constants 
for oxygen and nitrogen are so similar that the comparison would 
have very little bearing on the mixture equation of state. 

2. Masson and Dolley’ measured the compressibility of five 
mixtures of argon with ethylene, three mixtures of oxygen with 
ethylene, one of argon with oxygen, and of the three pure gases at 
24.95° C. and to 125 atmospheres. The gas compositions were ex- 
pressed as percentages by volume of the separate components each 
taken at 24.95° C. and 1 atmosphere. The volumes of each mixture 
were given in terms of the sum of the volumes of the component 
gases each measured at 24.95° C. and 1 atmosphere. 

3. Verschoyle** determined the compressibility of three mixtures 
of hydrogen and nitrogen and of the two pure gases at 0° and 20° C. 
to 200 atmospheres. The meaning of the term “percentage com- 
position” was not stated, so that it was assumed that this expression 
denoted percentages by volume of the separate components each 
taken at 20°C. For each mixture the unit of volume was the volume 
of that mixture at 0° C. and 1 atmosphere. 

4. Keyes and Burks’ *8 studied the compressibility of three 
mixtures of nitrogen and methane and of methane from 0° to 200° C. 
and to about 300 atmospheres. Smith and Taylor” using the same 
apparatus had previously studied the compressibility of pure nitrogen. 
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The composition of the mixtures were given in percentages by weight, 
and volumes in ccm. per gram. 

5. Bartlett, Cupples and Tremearne*® determined the compressi- 
bility of one mixture of hydrogen and nitrogen and of the pure gases 
from 0° to 400° C. and to 1000 atmospheres. These data are not 
treated in the present paper. 

6. Gibby, Tanner and Masson*® studied the compressibility of one 
mixture of helium and hydrogen and of the pure gases from 25 to 
175° C. and of several mixtures at 25° C. The original data are 
not given, but they list the parameters a and b of the relation: 


po= at bp 


for each isotherm of each mixture. The curve obtained by plotting 
the values of b against the composition of the corresponding mixture 
for 25° C. exhibits a maximum of about 80% hydrogen. Thus 
some mixtures of hydrogen and helium are less compressible than 
either of the pure gases. 

7. Scott measured the compressibility of three mixtures of hydro- 
gen and carbon monoxide and of the pure gases at 25° C. and to 170 
atmospheres. The meaning of “percentage composition” is not 
given but it is stated that the compositions were determined by 
chemical analyses. This may mean that the compositions are per- 
centages by volume at some low pressure or true molal percentages. 
In our treatment of these data it was assumed that compositions were 
expressed in mole percent, the difference between the two methods of 
giving compositions being probably within the experimental error of 
the measurements. Volumes were given in terms of the volumes of 
the mixtures at 25° C. and 1 atmosphere each of which was assumed 
to have the value 1.09155 under these conditions. 

In all of the measurements mentioned above temperatures were 
expressed in degrees centigrade and it is assumed that pressures were 
in normal atmospheres. 

In Table I is given a résumé of the methods of expressing composi- 
tions used by the four groups of investigators whose data are treated 
in the present paper, together with the corresponding mole fractions. 
The values of the molal volumes of the pure gases used for the data 
of Masson and Dolley and of Verschoyle were computed from the 
equation of state. 
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6. TREATMENT OF THE EXPERIMENTAL DatTa. 


The values of the constants of the Beattie-Bridgeman equation of 
state for several gases have been published*®* in units of atmospheres, 
liters per mole and degrees Kelvin. In Table II are given a résumé 
of the volume units in which the data of the four groups of observers 
are expressed and the conversion factor by which the density in 
moles per liter must be multiplied to reduce it to the units in which 
the data are given. 

In order to calculate the value of the conversion factor for Masson 
and Dolley’s measurements it is required to know the molal volumes 
of the pure gases at 24.95° C. and 1 atmosphere, and these values were 
computed from the equation of state for these gases. 

For the measurements of Verschoyle and of Scott it is necessary 
to know the molal volumes of the mixtures at 20° C and 1 atmosphere. 
and at 25° C. and 1 atmosphere respectively. These values were 
computed from the equation of state for mixtures. It is evident 
that this computed value depends on the method of combination of 
constants used. Hence in each case one factor was calculated using 
the linear method of equation (19), and one factor using the method 
of (19) for all constants except Bo, for which the Lorentz method (20) 
was used. The first factor was employed whenever linear combina- 
tion for Bo, was used as in Table IV and the second when Lorentz 
combination was used, Table V. 

The measurements of Keyes and Burks were expressed in volume 
units of ccm. per gram, and the calculations were made in the same 
units. 

In Table III are given the values of the equation of state constants 
for the pure gases and of the parameters for the mixtures calculated 
according to the relations (19). There is also given in each case the 
value of Bo, obtained by using Lorentz combination, equation (20). 
For the data of Keyes and Burks the weight fractions and the values 
of the constants in units of cem. per gram were substituted into (19). 
In the Lorentz combination of By for the latter data however the 
calculation was carried out using mole fractions and the values of 
By for the pure gases in cem. per mole. The value of Bo, so obtained 
was divided by the average molecular weight of the mixture given in 
Table I to obtain the value of Bow which can be used when volumes 
are in ccm. per gram. 

In Table IV are given the values of the virial coefficients Bn, yn 
and 6, of equation (18) for each of the pure gases and mixtures at 
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each of the temperatures for which compressibility data exist. The 
values of the parameters for the mixtures are those obtained by use 
of the method of combination (19) and appear in Table III. The 
virial coefficients are given first in units of atmospheres, liters per 
mole and degrees Kelvin and then in the volume units of the data. 
The conversion was made by dividing RT, Bn, yn and 6, by the 
appropriate factor given in Table II raised to the first, second, third 
and fourth powers respectively. 

In Table V the values of the virial coefficients are given when the 
method of combination (19) is used for all constants except Bo which 
is combined by the Lorentz method (20). The conversion factors 
used in this case for the data of Verschoyle and Scott are those listed 
in Table II under the heading ‘‘ Lorentz Bo;.”’ 

The virial coefficients were not computed for the data of Keyes 
and Burks since their measurements were made by the isometric 
method. The calculated values were computed accordingly, as 
described in an earlier paper.”® 

It is to be noted that there are no measurements of the compressi- 
bility of carbon monoxide that extend over a sufficiently wide tempera- 
ture and pressure range to permit the determination of the equation 
of state constants for this gas. Computations*® of the weight of a 
normal liter of carbon monoxide by use of the constants for nitrogen 
have confirmed the conclusion drawn from Langmuir’s* theory of 
isosters that the constants for carbon monoxide should not differ 
greatly from those for nitrogen. Accordingly, in the present com- 
parison the constants for nitrogen were used for carbon monoxide 
and this procedure is justified by the fact that to 170 atmospheres 
the deviation between the calculated and observed pressures for 
carbon monoxide are of the same order of magnitude as those for 
hydrogen (see Figure 21). 

A value of the pressure was calculated from the virial coefficients 
given in Table IV, and one from the values given in Table V for 
each observed point of Masson and Dolley, Verschoyle, and Scott. 
Corresponding computations were made with the aid of the values 
given in Table III for the results of Keyes and Burks. The per- 
centage deviations of the calculated from the observed pressures: 


: obs. pressure — caled. pressure 
% = 100 —> P 





obs. pressure 


were plotted against the observed pressures for each case. The 
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plots are given in Figures 1 to 24. Figures 1, 10, 14 and 21 exhibit 
the results for the pure gases. For each mixture there are two 
curves: one, marked “Lin,” representing the deviations when the 
method of combination (19) is used for all constants; and, the other, 
marked “Lor,” giving the deviation when the relations (19) are 
used for all of the constants except Bo which is combined according 
to the Lorentz method (20). On each mixture plot for the data of 
Masson and Dolley, Verschoyle, and Scott there also appears a 
curve, marked “‘Sum,” the ordinate of which at a given pressure is 
the sum (21%, + 22%:2) where 2; and 22 refer to the mole fraction of 
the Gases 1] and 2 in the mixture, and % and %,2 to the percentage 
deviations (from Figure 1, 10, 14 or 21) for the pure gases at that 
pressure. 

Although in some cases there were slight irregularities, smooth 
curves were drawn except for the data of Verschoyle. No curves 
are given for the argon-oxygen mixture since the three curves “ Lin,” 
“Lor,” and “Sum” practically coincide when plotted on the same 
scale as that used for the other figures. 
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7. GRAPHS OF THE DEVIATIONS OF THE PRESSURES CALCULATED BY 
THE BEATTIE-BRIDGEMAN EQUATION OF STATE FOR MIX- 
TURES FROM THE OBSERVED VALUES. 


Percentage deviation. 
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Pressure, atmospheres. 


FiGuRE 1. Masson and Dolley: Argon, oxygen and ethylene at 24.95° C. 
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Figure 2. Masson and Dolley: 75.17% argon, 24.83°7% ethylene at 24.95° C. 
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FIGURE 3. Masson and Dolley: 49.92% argon, 50.08% ethylene at 24.95° C. 
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FigurE 4. Masson and Dolley: 40.02% argon, 59.98% ethylene at 24.95° C. 
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Figure 5. Masson and Dolley: 29.18% argon, 70.82% ethylene at 24.95° C. 
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FiGuRE 6. Masson and Dolley: 9.89% argon, 90.11°% ethylene at 24.95° C. 
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FicurE 7. Masson and Dolley: 74.63% oxygen, 25.37; ethylene at 24.95°C, 
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FiGuRE 8. Masson and Dolley: 49.96% oxygen, 50.04% ethylene at 24.95° C. 
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FiGcurE 9. Masson and Dolley: 40.04% oxygen, 59.96% ethylene at 24.95° C. 
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Figure 10. Verschoyle: Hydrogen and nitrogen at 20° C. 
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FicgurE 11. Verschoyle: 74.98% hydrogen, 25.02% nitrogen at 20° C. 
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FiGgurRE 12. Verschoyle: 49.98°% hydrogen, 50.02% nitrogen at 20° C. 
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Figure 13. Verschoyle: 24.98% hydrogen, 75.02% nitrogen at 20° C. 
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Figure 14. Verschoyle: Hydrogen and nitrogen at O° C. 
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FIGURE 16. 
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Verschoyle: 74.98% hydrogen, 25.02°% nitrogen at 0° C. 
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FicgurE 17. Verschoyle: 24.98°% hydrogen, 75.02% nitrogen at 0° C. 
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Figure 18. Keyes and Burks: 56.66% nitrogen, 43.34% methane at 
0°, 100° and 200° C. 
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FicgurE 19. Keyes and Burks: 20.46% nitrogen, 79.54% methane at 
0°, 100° and 200° C. 


















































—0.5 





100 150 200 250 300 350 400 


Pressure, atmospheres. 


Keyes and Burke: 19.46% nitrogen, 80.54°> methane at 
0°, 100° and 200° C. 


Figure 20. 
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FiGuRE 21. Scott: Hydrogen and carbon monoxide at 25 C. 
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FiGuRE 22. Scott: 66.30, hydrogen, 33.7°% carbon monoxide at 25° C. 
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FIGURE 23. Scott: 51.7°> hydrogen, 48.3°% carbon monoxide at 25° C. 
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Figure 24, Scott: 33.1°% hydrogen, 66.9°% carbon monoxide at 25° C 
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8. DISCUSSION OF THE GRAPHS. 


We shall first consider the results exhibited in the graphs with a 
view to determining which of the methods of combination of the 
constants of the equation of state gives the best representation of the 
pressures of the mixtures. 

The constants a, b and c have in general so small an effect on the 
calculated pressure that the method of combination used for them is 
not important; hence throughout the present treatment linear com- 
bination is employed for each. Of the two methods of combination 
that could be used for Ao, linear or linear square root, the former is 
excluded because this would lead to zero or to negative intercepts 
for the expansion of gases on mixing at zero pressure; hence in all 
cases linear square root combination is used for Ao. For Bo three 
methods have been proposed: linear, linear square root, and Lorentz. 
The plots have been drawn to exhibit the results for these methods 
of combining Bo. 

For each mixture there is given two curves one representing the 
percentage deviation of the calculated from the observed pressure 
when linear combination is used for By and the other when Lorentz 
combination for Bo is used. In each case all of the other constants 
are combined according to the relations (19). The deviation curves 
resulting when linear square root combination is used for By would 
in each case be about as far above the “Lor” curve as the “Lor” 
curve lies above the “Lin” curve. In the first place we notice that 
the differences between the results obtained by using the linear or 
Lorentz method of combination for Bo are not great and that they 
increase with increasing differences in the values of By for the gases 
composing the mixture. 

The measurements of Keyes and Burke together with those of 
Smith and Taylor are the only set that extends over a sufficient 
temperature range to permit the determination of the values of the 
equation of state constants for the pure gases from data taken with 
the same apparatus as that employed for the mixtures. In Figures 
18, 19 and 20 are shown the deviations of the calculated from the 
observed pressures for three of the five temperatures, 50° and 150° 
not being included. The average deviations for all temperatures and 
for the entire pressure range are given in Table VI. There is a 
slight trend in favor of linear combination for Bo. 
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For the discussion of the results of Masson and Dolley, of Verschoyle, 
and of Scott from the standpoint of the methods of combination of 
constants we are not in so favorable a position since it is necessary 
to employ values of the constants for the various pure gases deter- 
mined from the data of other observers. It is generally true that 
two sets of compressibility data taken by different observers on the 
same gas do not agree with each other to the same order of accuracy 
as the consistency of each group of measurements within itself. For 
this reason the curve labeled “Sum” which is determined by >(2;%;) 
is added. It is not intended that this curve correct completely for 
the deviations of the equation of state from the compressibility data 
on the pure gases, but it is believed that unless otherwise noted 
weight should be given to the curve when discussing the deviation 
curves for the mixtures. For reasons given below the curves repre- 
senting the Masson and Dolley data are not referred to the “Sum” 
curve; but those for the results of Verschoyle and of Scott are treated 
from this standpoint. 

The values of the constants for ethylene determined® from Amagat’s 
data do not give a satisfactory representation of Masson and Dolley’s 
measurements on pure ethylene although they do lead to a good 
reproduction of the mixture data. This lack of agreement in the 
case of the pure ethylene may be due in part to the fact that the 
densities above 70 atmospheres at 25° C. are greater than the critical 
density and in part to the real disagreement between the data of 
Amagat and of Masson and Dolley. It will be noticed that the 
mixture equation gives a fairly good representation of the mixture 
containing as much as 90[% ethylene and this suggests that perhaps 
there may be an error in the pure ethylene data. In any case it is 
evident that very little weight can be given the “Sum” curve when 
discussing the mixture data and hence the curves are compared with 
the axis representing zero deviation. An examination of the curves 
given in Figures 2 to 9 from this standpoint shows a trend in favor of 
linear combination of Bo. 

The results of the treatment of the data of Verschoyle and of Scott 
are given in Figures 10 to 17 and Figures 21 to 24 respectively. In 
each case the deviation curves for the mixtures are to be compared 
with that labeled ‘“‘Sum.” There seems to be little choice in the 
method of combining By for these mixtures since the curve marked 
“Sum” lies about half way between the “Lin” and “ Lor” curves in 
every case except for the mixture containing 74.989, hydrogen— 
25.024, nitrogen when it coincides with the “Lor” curve. 
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The linear square root method of combination for Bo leads to a 
deviation curve lying above the “Lor” curve. But the deviation 
curves resulting when Lorentz combination is used lie almost entirely 
above the experimental results, and hence linear square root com- 
bination for Bo is not supported. 

It is possible that one method of combination will give the best 
results for one pair of substances and another method for another 
pair. If however we wish to use one method for all mixtures, the 
evidence is slightly in favor of linear combination for Bo; and in 
view of the much greater simplicity of this method of combination 
it is the one recommended. 

It seems reasonable to conclude that the Beattie-Bridgeman 
equation (18) for gas mixtures together with the method of combina- 
tion (19) gives a satisfactory over-all correlation of the data for the 
‘ases studied, though the agreement may not be within the experi- 
mental error in some cases. 

The question of whether or not the method of combination (19) 
reproduces the data on gas mixtures with sufficient accuracy for 
determining the general thermodynamic properties of mixtures 
cannot be treated in the present paper. Investigation by Lurie and 
Gillespie,'’* Gillespie,*® and Gillespie and Beattie“ indicate that the 
method of combination (19) does give satisfactory results in several 
cases. It is believed that the most severe test of an equation of 
state for mixtures is the representation of compressibility data. 
An equation that gives a fair reproduction of the pressure-volume- 
temperature data on gas mixtures will represent many other thermo- 
dynamic properties excellently. 


9, LINEAR COMBINATION OF CONSTANTS COMPARED WITH DALTON’S 
AND WITH AMAGAT'S Law. 


The method of combination (19) may be called linear combination 
of constants. It is of interest to determine whether this method 
leads to subsiantially better results than Dalton’s law of the additivity 
of pressures at constant volume and temperature or Amagat’s law 
of the additivity of volumes at constant pressure and temperature. 

Lurie and Gillespie,"? and Keyes and Burks have found that 
linear combination of constants in the Keves equation of state gives 
better results than Dalton’s or Amagat’s law for mixtures of nitrogen 
with ammonia and of nitrogen with methane respectively. Guillespie*® 
has shown that linear combination of constants in the Beattie- 
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Bridgeman equation of state gives a much better representation of 
the expansion of gases on mixing at low pressure than Amagat’s 
law for argon-ethylene, oxygen-ethylene, nitrogen-methane and 
hydrogen-nitrogen mixtures. 

Keyes and Burks do not discuss their compressibility data from the 
standpoint of Dalton’s law of additive pressures; they use the law 
of “partial pressure”’ in terms of weight fractions. They also give 
comparisons (their Figure 1) of the volumes calculated by use of 
Amagat’s law with the observed volumes for two of their mixtures. 
In Table VII are given the percentage deviations (read from their 
plot) of the volumes of the mixtures calculated by Amagat’s law 
from the observed volumes; and the percentage deviations of the 
pressures of the mixtures computed by linear combination of constants 
in the equation of state from the observed pressures, these latter 
deviations being determined from our plots 18 and 20. Linear 
combination of constants gives a better representation of the nitrogen- 
methane data than Amagat’s law; and, since Amagat’s law should 
give better results for these mixtures than Dalton’s law, it is safe to 
conclude that linear combination of constants gives a better repre- 
sentation than Dalton’s law. 

Masson and Dolley give for all of their mixtures the deviations 
resulting when Dalton’s law and when Amagat’s law are used. In 
Table VIII are given the average percentage deviations of the pres- 
sures of the mixtures calculated by Dalton’s law from the observed 
pressures, and of the volumes of the mixture calculated by Amagat’s 
law from the observed volumes. The deviations of the pressures 
calculated by linear combination of constants from the observed 
pressure are also given. The Dalton law and Amagat law percentages 
were obtained by dividing the deviations by the calculated pressures 
or volumes, whereas the linear combination percentage deviation were 
obtained by dividing by the observed pressures in the usual manner. 
Examination of Table VIII shows that linear combination of constants 
is far superior to either adding pressures or adding volumes. 

Neither Verschoyle nor Scott treated their results from the stand- 
point of Dalton’s or Amagat’s law. Calculation for several points 
shows that for these mixtures linear combination of constants gives 
the best results but the superiority of this method is not so pronounced 
as for those mixtures containing ethylene. 
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10. THE SECOND VIRIAL COEFFICIENT. 


Gibby, Tanner and Masson‘? have found a maximum in the second 
virial coefficient for helium-hydrogen mixtures at 25° C. Using the 
method of combination of constants (19) we can write the second 
virial coefficient 8, for a binary mixture in terms of the mole fractions: 


8, = RT|a,Bo. + (1 — 21) Boo] — [x1 Ao’? + (1 — 2) Agi? 
R. | 
a pe ees +(1—21)c.] (22) 


Differentiating and equating to zero we find that 8B, has a maximum 
or minimum value at 





bs 2 2(Ao,!/? — Age! 3) si 
where 8; and #2 are the virial coefficients for the pure Gases 1 and 2 
respectively. Moreover, this value of x; must correspond to a 
maximum value of 8, since the curvature 


() = —? (A9,'/2 ‘aa Age/?)? 24) 
T 





dx? 


is always negative whatever the values of Ao; and Age. 

The existing experimental data on binary mixtures indicate that 
the B, vs. 2, curve does not exhibit a maximum, except perhaps in 
certain temperature regions. From equation (23) we see that no 
maximum will exist for values of 2; between zero and unity unless 
the numerical values of the last term of (23) is less than 14. Sub- 
stitution of the numerical values of the constants for hydrogen and 
helium into (23) gives a calculated maximum at about 50% hydrogen. 
The plot given by Gibby, Tanner and Masson indicates a maximum 
for B, at about 80% hydrogen. In a later publication,® Tanner and 
Masson “restate the data” of their earlier paper and give a corrected 
plot which indicates a maximum at about 70° hydrogen. The 
curvature is so slight that the maximum cannot be located with 
exactness. The lack of quantitative agreement between the calculated 
and the observed values of the maximum is probably due to one or 
both of two causes: (1) Deviation of the equation of state of hydrogen 
or helium from the data of Gibby, Tanner and Masson on the pure 
gases, and (2) the fact that these investigators used only two terms 
of the virial to represent their isotherms, whereas the equation of 
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state virial contains four terms for hydrogen and three for helium. 

Equation (23) calculates a maximum in the 6, curve for hydrogen- 
nitrogen mixtures at 97°7 hydrogen for 20° C. but no maximum at 
0° C. Verschoyle’s curves were drawn without a maximum, but the 
20° curve almost admits of a maximum. 

Gibby, Tanner and Masson, and Verschoyle each expanded pV in an 
ascending power series of the pressure but the second coefficient of 
this expansion should correspond to the 8, of equation (22) divided 
by RT. 

11. SUMMARY. 

The Beattie-Bridgeman equation of state (18) has been applied to 
the calculation of the pressures of gaseous mixtures of argon with 
ethylene, oxygen with ethylene, hydrogen with nitrogen, nitrogen 
with methane and hydrogen with carbon monoxide. 

The equation of state constants a, b and c, which do not, in general, 
have a great effect on the calculated pressure were combined linearly 
with the composition of the mixture. In accordance with the results 
of several previous investigations, lmear square root combination 
was used for Ap. The relative merits of linear, Lorentz, and linear 
square root combination for By were studied by the use of deviation 
curves for all of the mixtures. On the whole, linear combination 
for By was slightly favored over the Lorentz method, while the linear 
square root relation gave the least satisfactory results. It is pointed 
out that there is little theoretical evidence in favor of the Lorentz 
method. 

Linear combination for each of the constants R, A,o'’*, a, Bo, b and 
c of the Beattie-Bridgeman equation of state gives a satisfactory over- 
all correlation of the compressibility data on gas mixtures, although 
in some cases the representation of the data may not be within the 
experimental error. 

For all of the mixtures studied, linear combination of the constants 
of the equation of state gives a much better representation of the 
experimental data than either Dalton’s law or Amagat’s law. 

It had been experimentally found that certain mixtures of helium 
and hydrogen have a larger second virial coefficient than the pure 
gases. ‘The equation of state confirms the existence of such a maxi- 
mum and predicts a similar phenomenon for other mixtures within 
certain temperature ranges. 

The values of the equation of state constants for nitrogen give a 
satisfactory representation of the compressibility of carbon monoxide 
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to 170 atmospheres at 25° C., a result which is not unexpected in 
view of Langmuir’s theory of isosters. 
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